Necessary conditions for accurate computations of three-body partial decay widths 
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The partial width for decay of a resonance into three fragments is largely determined at distances 
where the energy is smaller than the effective potential producing the corresponding wave function. 
At short distances the many-body properties are accounted for by preformation or spectroscopic 
factors. We use the adiabatic expansion method combined with the WKB approximation to obtain 
the indispensable cluster model wave functions at intermediate and larger distances. We test the 
concept by deriving conditions for the minimal basis expressed in terms of partial waves and radial 
nodes. We compare results for different effective interactions and methods. Agreement is found with 
experimental values for a sufficiently large basis. We illustrate the ideas with realistic examples from 
a-emission of 12 C and two-proton emission of 17 Ne. Basis requirements for accurate momentum 
distributions are briefly discussed. 

PACS numbers: 21.45.-v, 31.15.xj, 25.70.Ef 



I. INTRODUCTION 

The spectrum of a given many-body quantum system 
provides a set of characteristic observables: the energies 
and the widths. All quantum states, except perhaps the 
ground state, decay if sufficient time is available. For 
bound states, where the total energy is less than all 
thresholds for division into subsystems, only electromag- 
netic and particle transforming decays are possible. The 
first of these decays maintains the identity of the con- 
stituent particles in contrast to the latter exemplified by 
beta-decay where neutrons are transformed into protons 
or vice versa. However, if the energy is sufficiently high, 
the system can also decay into subsystems while main- 
taining the identity of the constituent particles. In this 
work, we shall concentrate on decays where the energy 
allows such fragmentation of the initial system. 

In general we then have a many-body continuum prob- 
lem. For nuclei the simplest final state consists of two 
fragments, e.g. two fission fragments, or a daughter nu- 
cleus plus a nucleon or an a-particle [1]. If both ini- 
tial and final states are completely specified, energy and 
momentum conservation determine the relative kinetic 
energy between the two outgoing particles. The decay 
rate (or the width found by multiplying the rate by K) 
is obtained from preformation, or spectroscopic factors, 
combined with the probability for tunneling through the 
barrier in the relative potential created by the two-body 
interaction. This barrier separates the short-distance ini- 
tial many-body state from the large-distance final two- 
body state [2]. The many-body problem is reduced to 
a two-body problem where only the two particles found 
after the decay appear. 

Larger widths may be found in more elaborate mod- 
els exploiting different, perhaps virtual, configurations 
resulting in a coupled channels problem [3]. Such rel- 
atively simple two-body decays have been studied from 
the beginning of the history of quantum mechanics, ex- 
emplified by a-emission [1]. Similar processes vary from 



statistical emission of nucleons above the nucleon sepa- 
ration energy for ordinary nuclei [4], to (almost) instant 
decay outside the neutron dripline and to proton and a- 
emission outside the proton dripline [5, 6]. 

In this work we shall consider decay processes where 
three fragments are found in the final state. This is the 
simplest, yet not understood, extension of the concept of 
two-body decay [7]. Furthermore, to limit the number 
of possible final states we assume that the three-body 
threshold is lower than any other threshold. Energy and 
momentum conservation still provide constraints, but the 
internal distribution of the total momentum and energy 
between the three fragments is not decided by these 
conservation laws. These momentum and energy dis- 
tributions are observables carrying detailed information 
about the initial state and the process. Reliable com- 
putations require accurate determination of the large- 
distance properties of the cluster wave functions [8, 9]. 
On the other hand, the decay width is an average quan- 
tity which is very sensitive to properties of the potential 
barrier, but in analogy to two-body decay, it is deter- 
mined by the effective barrier at small and intermediate 
distances. 

Nowadays, different methods to compute the partial 
decay width into three specified fragments are already 
available [9-11]. However, the conditions for their relia- 
bility and suitability are not well established, and each 
method is most often rather tested on the individual sys- 
tems under investigation. The discussions comparing the 
different methods are deceivingly mixing effects of choices 
of (i) degrees of freedom, (ii) interactions, (hi) theoreti- 
cal method, and (iv) numerical convergence. Untangling 
these effects is badly needed to formulate necessary con- 
ditions for accurate computations. Benchmark compu- 
tations for precisely specified systems and interactions 
would be valuable as test criteria for reliability of the 
methods. 

The purpose of the present paper is to provide veri- 
fiable simple but revealing test examples. To do this it 



2 



is necessary to separate and assess the impact of each 
of the effects (i)-(iv). We first explain why the basic in- 
gredient for partial three-body decay widths necessarily 
must be a three-body cluster model. We shall formu- 
late necessary conditions for accurate three-body com- 
putations, and document by numerical applications on 
realistic three-body decaying systems. It is crucial that 
extensions to include more complicated effects are built 
on methods which are established as accurate. 

The basic concepts for three-body decay widths are in 
section 2 described and tested against measured results. 
In section 3 we give analytical estimates of the quantities 
characterizing the crucial potentials and the basis size 
needed for accurate computations. In section 4 we test 
the estimates with the clean example of the Hoyle reso- 
nance in 12 C. We then discuss in section 5 how the com- 
puted widths depend on the basic interactions and the 
available methods. We also briefly discuss the more se- 
vere accuracy requirements for momentum distributions 
determined at larger distances. Finally section 6 contains 
summary and conclusions. 



II. FRAMEWORK 

Different definitions exist of the decay width, e.g. in 
terms of cross sections [12], or phase shifts and S-matrix 
poles [13], or eigenvalues of the Hamiltonian H [14]. We 
shall use the width of a given resonance defined as minus 
two times the imaginary part of the generalized eigen- 
value of the Hamiltonian. This is equivalent to a pole 
in the 5-matrix at the complex momentum correspond- 
ing precisely to that eigenvalue. This definition is only 
complete after specification of the degrees of freedom con- 
tained in H. The full many-body Hamiltonian would give 
the total width, while confined to three particles in the 
final state the result is an approximation to the corre- 
sponding partial decay width. 



A. Basic concept 

The important issue in three-body decays is that the 
many-body degrees of freedom must be (re)organized into 
intrinsic and relative cluster coordinates. This division is 
usually not meaningful at small distances when all par- 
ticles are close and within the nuclear volume. In con- 
trast, this is the only meaningful division at large dis- 
tances where three free particles are present. The quan- 
tum mechanical wave function must reflect this transition 
and the three-body structure at large distance must un- 
fold into a many-body structure at small distance. To 
get a correct partial three-body width, effects from the 
small distances must be incorporated, e.g. by a preforma- 
tion factor and the assumption of an artificial attractive 
pocket designed to provide both the correct energy and 
the resonance small-distance boundary condition. This 



treatment is precisely as for a-emission in the classical 
Gamow theory [1, 2, 7]. 

In practice the many-body problem is therefore trans- 
formed into a three-body problem for all distances. The 
potential pocket in three-body coordinates has to be 
added by hand, unless the chosen two-body interactions 
already are sufficient. It is crucial to have the correct 
three-body resonance energy as evident from the expo- 
nential energy dependence of the width determined from 
tunneling probability through any barrier. Thus we have 
to insist on a practical method to adjust the energy to 
the correct value, e.g. by use of a short-range three-body 
interaction [15]. This separates the model dependence of 
short-distance many-body structure from effects of the 
three-body cluster model at intermediate distances. In 
Gamow theory this is achieved by adding an attractive a- 
daughter interaction, e.g. a square-well or Woods-Saxon 
potential with a radius about the size of the nucleus. 

All methods to compute partial three-body decay 
widths must at some point address this separation of 
distances, degrees of freedom and related effective inter- 
actions. However, the procedures to reach the reduc- 
tion into the three-body structure are rather different. 
Microscopic derivations start from nucleon-nucleon in- 
teractions and integrate away the unwanted degrees of 
freedom while simultaneously leaving corresponding ef- 
fective interactions [16-20]. The mean-field model is also 
used for two-proton decay in the formulation named the 
Gamow shell-model where the nucleon degrees of free- 
dom are present and the interactions are phenomenolog- 
ical adjusted within the mean-field [21-23]. 

If the three-body cluster model is assumed from the 
beginning for all distances, the corresponding two-body 
interactions are usually obtained from phenomenological 
adjustments to two-body data. No matter how the effec- 
tive two-body interactions are found, they are much more 
important in three-body than in two-body decays. For 
example for a-cmission, the distance between a-particle 
and daughter quickly leaves only non-vanishing contri- 
butions from Coulomb and centrifugal forces [8, 9]. For 
three-body decay, where two particles stay close while the 
third particle moves away, the short-range interactions 
contribute much more to the properties of the decisive 
barrier [10]. 

The crux of the matter is then to deduce the partial de- 
cay width from the three-body cluster model for a given 
resonance energy. Here several methods have been em- 
ployed. The practical and experiment oriented method 
is to use elastic scattering cross sections as function of 
energy [12, 13]. A peak is then related to a resonance 
and its width is the resonance width. However, this is 
not practical for collisions of more than two particles, 
and furthermore uncertainties arise from corrections due 
to phase space distortion, broad peaks, overlapping res- 
onances, background contributions, etc. We prefer the 
more mathematical definitions of a complex pole in the 
S- matrix [13], or the energy derivative of the scattering 
phase shift while crossing through tt/2 [12, 13], or equiv- 
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alently minus twice the imaginary part of the eigenvalue 
of the hamiltonian [14], or equivalently related to the 
solution of complex energy with an outgoing flux in all 
channels. 

The numerical method to compute the complex reso- 
nance energy has to allow for complex energy solutions 
and for example insist on only an outgoing flux in all 
channels [24]. Equivalently all coordinates can be ro- 
tated a given angle into the complex plane which turn 
the outgoing flux solution into an exponentially falling 
solution at large distances precisely as for bound states 
[25, 26]. Another method exploits analytic continuity of 
the interactions by varying with respect to for example 
the strength parameter [27, 28]. These methods are es- 
sentially all equivalent and produce identical results for 
identical interactions and the same Hilbert space. 

A semiclassical perturbative method has been abun- 
dantly employed recently [11, 29, 30]. It assumes that 
the width is small and related to the outgoing flux for the 
solution to the Sehrodinger equation obtained by confine- 
ment to a box of finite extension. This method is pre- 
sumably also equivalent to the other methods for narrow 
resonances when the perturbation assumption is valid. 



B. Testing the concept 

The bare problem can be envisaged as a potential and 
a corresponding wave function with resonance character 
conveniently described by a complex energy and only an 
outgoing flux. Real (En) and imaginary (£/) parts of 
the energy are closely linked through the potential, e.g. 
adding a small positive potential at small distance would 
increase En, leave the potential barrier unchanged, and 
increase \Ej\ exponentially. We shall test the idea with 
the simplest computation of the width, T = —2Ej, i.e. 
by use of the WKB-approximation. 

However, to do this, a potential defined as function of 
a generalized radial coordinate is required. The coordi- 
nates for three particles then should be combined into one 
overall important coordinate p which on its own should 
be able to describe the process. We want to maintain the 
intuitive understanding that small p means small physi- 
cal distances between all three particle at the same time. 
Vice versa, large p should be able to describe the struc- 
ture after fragmentation into three particles. To give a 
precise and physically meaningful definition of p in terms 
of the particle coordinates we assume a quadratic depen- 
dence on interparticle distances. Then the hyperradius p 
is unique apart from the choice of mass weighting, where 
we use mMp 2 — X)i<j m i m j(^i — ?j) ■> with m, and r*j as 
mass and coordinate of particle i, and M = mi+ra 2 +ni3 
[31, 32]. In our calculations we take the normalization 
mass to equal to the nucleon mass. 

We shall use p as the generalized coordinate. In the 
spirit of the generator coordinate method, where the en- 
ergy is calculated as function of one parameter, we solve 
the Sehrodinger or Faddeev equations for each value of p 



between zero and p ma x- The hyperangular degree of free- 
dom is automatically quantized as the solution to the cor- 
responding Sehrodinger equation, which for each p takes 
the form: 



A 2 $ n ( P ,n) + 



2mp 2 



v(x)^„( P ,n) = x n (p)^ n (p,n), (i) 



where A is the hyperangular operator [32] , V is the sum 
of the three two-body potentials, and f2 represents the 
five hyperangles. 

If we now use the complete set of solutions for each p as 
basis (* = E„/n(p)$n(p,fi)) we arrive at the hy- 
perspherical adiabatic expansion method, which reduces 
the problem to the following set of coupled differential 
equations in p: 

cP_ 

dp 2 n? p* 

■ X ( - ' ' . '< ' ~r ■ ! /'»'(/') ;:(s < 2 > 



?T77 1 

+ ^(V 3b (p)-E) + -[X n (p) + 
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Up) 



dp 



where E is the three-body energy, V36 is a three-body 
potential used for fine-tuning, and the functions P nn > and 
Qnn' are given for instance in [32]. 

These equations contain the effective adiabatic poten- 
tials, which take the form: 



VeS(p) 



h 2 A„Qo) + 15/4 
2m p 2 



V 3b (p). 



(3) 



The number of differential equations (or of adiabatic po- 
tentials) is equal to the size of the basis for each p. This 
basis is unique in the description of multifragmentation, 
because it is the only representation that maps fragmen- 
tation theory onto a set of coupled channel differential 
equations from two-body reaction theory [31, 33]. 

To fix the concept we shall first only use the lowest of 
these adiabatic potentials which depend on angular mo- 
mentum and parity of the three-body system. The three- 
body configurations change in a non-trivial manner from 
small to large hyperradii. The total energy, apart from 
the kinetic energy related to variation of p, is minimized 
for each p. However, p does not uniquely determine even 
the geometric configuration. The same p is related to 
continuously differing combinations of distances between 
the particles, e.g. a small distance between two particles 
and a large distance to the third particle, or equal dis- 
tance between all particles, etc. Thus, it is not a priori 
obvious why this should be a good choice of coordinates 
for an efficient description of these decays. Other config- 
urations might be important but this would be reflected 
in a finite population of the complete set of higher-lying 
adiabatic potentials. 

In Fig.l we show the lowest adiabatic potentials 
for various angular momenta and parities J 71 ' in 12 C 
(a+a+a). The behavior of the three body systems is 
unpredictable from J 71 " alone. A minimum at small dis- 
tances indicates a substantial amount of cluster struc- 
ture and vice versa, no minimum indicates dominance 
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FIG. 1: The computed lowest adiabatic potential as function 
of hyperradius for a number of resonances of 12 C (a+a+a). 
The horizontal lines mark the resonance energies measured 
above the three-body threshold. The corresponding excita- 
tion energies are given above each of the panels [34]. 



of many-body non-cluster structure often referred to as 
shell model structure. In all cases, after decay, the three 
alpha particles must emerge outside the barrier. The 
physical meaning of this description is different from the 
intuitive perception of a mean-field shell model where a 
monotonous dependence on excitation energy and angu- 
lar momentum would be expected. Here the three-body 
structure and the properties of the two-body interactions 
arc crucial and capable of changing the expected order- 
ing. 

The partial decay width must sensitively depend on 
the energy and the properties of the barrier which must 
be crossed. The complex energy solutions of the hyper- 
radial equation give the widths, which we here estimate 
by the WKB tunneling probabilities through the one- 
dimensional potential barriers. We fix the real part of 
the energies equal to the measured values whenever they 
are available. This can be achieved by using an appropri- 
ate three-body short-range attractive potential (V-^p) in 
Eq.(2)). A simple form of this is employed in Fig.l, where 
we take a square well potential of radius 4 fm and a depth 
adjusted to give the experimental resonance energies. 

The WKB estimates are obtained from the action inte- 
gral between the classical turning points determined by 
the real part of the energy (E). The transmission coeffi- 
cient is then given by: 
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FIG. 2: The computed widths compared to the measured val- 
ues for different resonances of 12 C with excitation energy E* 
[34] . The WKB approximation is used for tunneling through 
the lowest adiabatic potential where the inner turning point 
is fixed to a hyperradius of p = 4 fm. 



where pi and p t are the inner and outer classical turning 
points defining the distance through the barrier. Once 
the transmission coefficient is computed, the decay con- 
stant can be obtained as T/h=fT, where / is the knock- 
ing rate. 

The results are shown in Fig. 2 where only a factor of 
2-5 remains to match the measured values which vary 
by about 5 orders of magnitude with only a variation 
of 7 MeV in excitation energy E* . The smallest widths 
occur for the smallest angular momenta but the largest 
widths do not occur for the largest angular momenta. 
The only exception is the lowest 1 + state where the com- 
puted width is too small by about a factor of 100. The 
reason is that the large-distance tail of the correspond- 
ing potential sensitively contributes to the width for this 
state, see Fig.l. The basis size should then in this case 
have been larger than used to obtain this estimate. 

In more accurate computations, as discussed in details 
in this report, the inner turning point would also be larger 
implying larger computed widths. This would in turn in- 
dicate reduction due to preformation, or spectroscopic, 
factors, arising from a short-distance structure deviating 
from that of three a-particles. In any case, we can con- 
clude that the use of the hyperradius as the generating 
coordinate is valid to account for the main dependencies 
of the resonance widths on the excitation energy and an- 
gular momentum and parity of the system. 



III. CRUCIAL INGREDIENTS 

The exponential dependence of the widths on the prop- 
erties of the confining barrier immediately emphasizes 
that the most crucial ingredient is the potential barrier 
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which must be accurately determined between classical 
turning points. Any inaccuracies would be exponentially 
enhanced in the width computation. Thus, it is essential 
first to know the coordinate region where accuracy is in- 
dispensable which means that the turning points should 
be found. Second, it is crucial to control the numerical 
technique responsible for the computation of the poten- 
tial barriers in this coordinate range. 



A. Turning point estimates 

The inner turning point is at a distance where the 
short-range interaction is most important. It is sensi- 
tive to the details of the potential and the energy and 
quantum numbers of the resonance. However, the short 
distance allows accurate computation without too many 
difficulties, provided that the cluster division is assumed 
and the many-body problem reduced to that of three 
particles. 

Then the outer classical turning point is an estimate 
of the largest distance needed in accurate computations. 
If we assume that only repulsive Coulomb interactions 
remain at these distances we can use hyperradial coordi- 
nates and find the turning point p t with the formalism 
developed in [35]. 

When a direct decay of the three-body system is as- 
sumed, all relative two-body distances scale proportion- 
ally. We can then minimize the action integral and arrive 
at the general expression [35] 

2 / \ 3/2 

5>^ fe ) 2/3 K™fc) 1/3 (5) 

\i<k ) 



Pt = 



EVmM 



where ezi is the charge of particle i and E is the energy 
of the three-body resonance which at the turning point 
equals the potential energy. 

The partial waves necessary to describe this distance 
can be estimated from the fact that at the turning point 
the value of the potentials can not be higher than the 
three-body energy E. Therefore, the centrifugal barrier 
must be smaller than E at the turning point, i.e. 



ft 2 , 4(4 + 1) + iy{iy + l) < E ^ 

2m x 2 y 2 



(6) 



where x and y are the usual Jacobi coordinates, and £ x 
and £ y are relative angular momenta related respectively 
to the distance between two particles and their center of 
mass and the third particle. Assuming that both the to- 
tal angular momentum and the intrinsic spins are small, 



we have 



„ and for sufficiently large values of 



and £ y we can take 4(4 + 1) w £ 2 X and £ y (£ y + 1) w £ y . 
Furthermore x is proportional to the distance between 
two of the particles, say 1 and 2, and for a direct decay 
all distances between pairs of particles are similar. This 
leads to 
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with [i\2 = mi+m 2 ' f rom which one gets 



i < ?! 



^2EmMm 1 m 2 m d , . (8) 



Two limits of Eqs.(5) and (8) are useful in practice. 
When all particles have the same masses mo and charges 
z we find 
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(9) 



When particle 3 has mass 777,3 and charge z 3 much larger 
than for the other two particles, i.e. mi w m 2 -C 777,3 an d 
Z\ w z 2 <C z 3 , we obtain 
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z\z z e' 
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8m 1 
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£ y < 2ziz 3 e" 



mi 
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(10) 



Notice here that the estimates of lengths in hypersphcr- 
ical coordinates always are combined with the normal- 
ization mass m, e.g. only pt\fm is expressed in terms 
of physical parameters like particle energies, masses and 
charges. In contrast dimensionless quantities like the an- 
gular momenta in Eqs.(9) and (10) do not involve the 
arbitrary mass m. 

Instead of following a path where all distances 
scale proportionally it could be advantageous to tunnel 
through the barrier by exploiting the two-body attrac- 
tion between two particles while the third particle moves 
away. This is a sequential decay where the geometry re- 
duces the Coulomb interactions and supplies additional 
energy from the short-range interaction. In total the bar- 
rier could be substantially smaller both in height and 
width. If the intermediate two-body configuration car- 
ries an energy Ei 2 we find 



Pt 



(E — £12) 
y < z z {z x + z 2 )e 



Z3,{zi+z 2 )e 2 m 3 (mi+m 2 ) 



mM 



2 / 27773(7771 + m 2 ) 
Mh 2 (E-E 12 ) ■ 



(11) 
(12) 



Where we have used that for sequential decay large values 
of p imply p rts y, and the relative distances r 3i and r 32 
are similar to n 2 ,3 = y\/mj p 12 , 3 w p^m/p,i 2 ,3- The 
angular momentum estimate is an upper limit because 
the two-body attraction can typically only be exploited 
for one or very few given small values of £ x which by 
angular momentum conservation also forces £ y to have a 
small value. 

In the two limits of the same masses (mo) and charges 
(zo), and one mass and charge (7773 and 03) much larger 
than the other two (mi « to 2 and z\ ~ z 2 ), we find that 
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(13) 
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The balance between the smaller Coulomb energy and 
the additional two-body energy _E 12 determines whether 
the tunneling process is direct (all pairs of particles are 
outside their short-range attraction at the turning point 
and all distances scale proportionally) or sequential (the 
attraction from one or more pairs of particles is used as 
vehicle to speed up the process by tunneling through a 
smaller barrier). Different paths can contribute to the 
same decay process. The dominating path would typi- 
cally correspond to the smallest outer turning point p t - 

B. Basis requirements 

The width can only be accurately computed if the po- 
tential and the related wave functions are accurate up 
to at least the true outer turning point. For this it is 
important that the Hilbert space is sufficient to allow 
the system to choose the optimum path through an ac- 
curately determined barrier. This in turn can only be 
achieved when the numerical procedure allows a precise 
coherent description of all two and three-body interme- 
diate configurations. 

For coherently contributing two-body substructures 
this is only ensured by Faddeev, or Faddeev-like, decom- 
positions. One set of Jacobi coordinates is obviously sim- 
pler than including all three Faddeev components but it 
also has a tremendous disadvantage for systems where 
more than one two-body subsystem simultaneously have 
bound or nearly bound states, and especially when rela- 
tively large distances are important. The Efimov effect 
is for example completely excluded from the description 
in a basis of one Jacobi set of coordinates. 

In many computations basis expansions arc exploited. 
With the hyperradius as coordinate the natural basis is 
hyperharmonics for each Jacobi set in the basis. To de- 
scribe structures varying over distances comparable to 
the range of the two-body interaction, b\2, the basis 
must on average contain a few points within that dis- 
ta nce. F or a given p the two-body distance is r\i = 
\Jm/ ' p,\ix (x = psina), which means that a proper de- 
scription of the internal two-body structures for large p's 
requires ba sis term s with nodes within the small region 
< a < \J \i\ilm byil p. The number of nodes between 
an 7r/2 of a given hypcrspherical harmonic is given by 
the principal quantum number n of the Jacobi polyno- 
mial contained in it. This means that a typical separation 
between nodes is 7r/(2n). Therefore, to have at least one 
node below \fpY2Jm b\2/ P one needs 

7r p I m 

n>~T-\ ; ( 15 ) 

2612VM12 

which implies that the maximum hypcrmomentum quan- 
tum number K max (= 2n+£ x +£ y ) should at least exceed 
a lower limit, i.e. 

K max > ^-<h^ + t x +ty , (16) 
bi2 V M12 



where we assumed given hyperradius and partial angu- 
lar momenta. Thus the number of basis functions, or 
K max , must increase proportional to p. The proportion- 
ality factor providing the unit of p for this estimate is 
inversely proportional to the range of the two-body in- 
teraction which initially is responsible for the structures 
we try to describe. The minimum basis size can only be 
found in practice by numerical calculations. 

When the turning point, p — pt, is needed we can 
insert the expression in Eq.(5) in Eq.(16). Then only 
the physical parameters remain in the estimate of K max . 
Furthermore, to relate pt to the physical size of the sys- 
tem the definition of the hyperradius can be used, which 
leads to mp\ w r1 v ^2 i<k m t m k /M where r\ v is an ap- 
propriate average of the distance, (r,— r^) 2 , between pairs 
of particles. 



IV. A TEST CASE: THE HOYLE STATE IN 12 C. 

We consider here the Hoyle state, i.e. the first + res- 
onance in 12 C described in terms of three a-particles, 
which is the only possible particle decay mode. The in- 
trinsic a-particle spins are zero and the wave function 
must be symmetric with respect to all interchanges of 
pairs of particles. This limits the number of necessary 
partial waves and the -ftT ma:r -dependence of the effective 
potential and wave function can conveniently be investi- 
gated. 

The well-studied converged result is known to be a nar- 
row resonance with a width of only 8.5 eV at an energy of 
0.38 MeV. A detailed calculation using the hypersphcri- 
cal adiabatic expansion method can be found in [9] where 
the structure is shown to arise essentially from only one 
of the adiabatic potentials. 

Following the estimates in Eq. (9) where a direct decay 
is assumed, we obtain an estimated value for the outer 
turning point of pt ~ 90 fm, while the maximum val- 
ues for the two-body relative orbital angular momenta 
are £ x w £ y w 7. However, the proportional scaling of all 
distances is inefficient because the two-body attraction is 
not exploited while the full Coulomb repulsion is encoun- 
tered. By using Eq.(13) instead we find p t ~ 67 fm with 
E12 = 0.1 MeV as the energy of the 8 Be subsystem, and 
£ y i=a 8. The coherence of all three subsystems reduces 
the turning point and keeps the number of components. 
Then the barrier is reduced and this escape mechanism is 
preferred. For the second case (p t ~ 67 fm and £ y w 8), 
and following (16), we get that K max has to be about 80 
in order to ensure a proper treatment of the 8 Be structure 
up to the outer turning point. 

In the left part of Fig. 3, we show the lowest adiabatic 
effective potential for the + states in 12 C for different 
values of K max . The dashed straight line indicates the 
energy of the resonance. In this calculation all the pos- 
sible values of l x and of £ y consistent with a given K max 
have been included. The curve quoted as full has been 
computed with K max —150 and values of l x and l y up 
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FIG. 3: The dominating effective adiabatic potential for the 
lowest + resonance in 12 C (a+a+a) for (a) different values 
of K m ax and (b) different values of partial waves in the ex- 
pansion in terms of HH. The dashed straight line indicates 
the energy of the resonance. In (a) all the the possible val- 
ues of £ x and £ y consistent with K max have been included. 
The curve called "full" includes £ x = £ y up to 12 only and 
K max — 150, but with K max increased up to 500 for some of 
the components. In (b) K max is taken equal to 150. 



to 12 only, but increasing K max up to 500 for the most 
contributing components. As seen in the figure, the com- 
puted potential barrier (and the outer turning point) 
changes dramatically with K maxi and, in agreement with 
the estimate of Eq.(16), a K max value of at least 80 is 
needed in order to obtained a converged potential up to 
the outer turning point, which is found to be 60 fm (also 
in agreement with the estimated value). 

In the right part of the figure we show the same effec- 
tive potential for a fixed value of K max (we have taken 
150), but where the values of i x and £ y are progressively 
increased. It is clear from the figure that too small val- 
ues of l x and l v produce a too small potential barrier. 
As estimated above, values of i x and £ y of at least 8 are 
needed to match the thick curve, which is the one plotted 
in the right part with K max = 150 and with values of £ x 
and l y of up to 74. 

It is then clear from Fig. 3 that simultaneous use of ap- 
propriate K max and maximum values of £ x and £ y are suf- 
ficient to reproduce accurately the potential barrier up to 
the outer turning point. This is the decisive region of the 
potential determining the width of a given resonance. It 
is important to note that while a too small value of K max 
overestimates the potential barrier, too small values of £ x 
and £ y underestimate it. Therefore both effects tend to 
compensate each other, in such a way that a poor calcu- 
lation using too low K max and too few £ x and £ y could 
however luckily give a computed width not too far from 
the correct one. 

The dependence of the width on the basis size is illumi- 
nating. In table I we give the computed WKB widths for 



TABLE I: WKB estimates for the width of the lowest + 
resonance in 12 C for the seven effective potentials shown in 
Fig. 3a (left part of the table) and for the seven potentials 
shown in Fig. 3b (right part of the table). The column labeled 
as rj, ert refers to the widths obtained with the perturbative 
method described in [11]. All the widths are given in MeV. 
The experimental value is 8.5 ■ 10 -6 MeV. 
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the different effective potentials shown in Fig. 3a (left part 
of the table) and Fig. 3b (right part of the table). The 
conclusion is striking, an insufficient basis (small K max 
or small maximum values of £ x and £ y ) easily leads to 
widths overestimated or underestimated by 3 — 5 orders 
of magnitudes. However the converged result agrees rea- 
sonably well with the experimental value of 8.5 eV. One 
has to keep in mind that only the uncertainty arising 
from the different possible choices of the knocking rate 
when computing the WKB width can easily produce vari- 
ations in the width of up to a factor of 2 or 3. In our 
calculations we have taken the knocking rate equal to 
the energy of the resonance divided by the Plank con- 
stant (E = hv). Furthermore, in this estimate we have 
neglected the higher lying adiabatic potentials which typ- 
ically contribute about 10-15% of the wave function. In 
any case, the width is exponentially depending on the 
barrier, and we have therefore demonstrated how catas- 
trophic it is to use an insufficient basis. 

A sufficiently large basis for three-body quantities is 
already important for structure computations of weakly 
bound halo systems [36]. Excitations into the contin- 
uum of two-neutron halos, where the Coulomb interac- 
tion is absent, give more sensitivity and the opportunity 
to compare computations with measured strength func- 
tions. The inadequacy of a small basis was exhibited in 
[37] for the 1" excitation of the 6 He ground state. This 
was further emphasized in the fairly successful prediction 
of the similar 1~ strength function for 11 Li in [38] which 
later on was accurately measured in [39]. 

V. INTERACTIONS AND METHODS 

For a-emission it is well-known that the long-range in- 
teractions (as the Coulomb potential) and the centrifu- 
gal forces, by far have the largest influence on the decay 
rates. This is less obvious for three-body decays since the 
potential barriers depend on the two-body interactions. 
On the other hand the decay rate should be independent 
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FIG. 4: Lowest adiabatic effective potential for the 3/2 _ reso- 
nance in 17 Ne. Except for the last curve (thick-dot), the p- 15 
interaction is the one described in [10]. The different curves 
correspond to calculations with different values of K max and 
relative two-body angular momenta. The thin-dotted curve 
corresponds to the calculation presented in [10]. In the last 
(thick-dotted) curve, the p- ls O potential and K max value of 
rcf.[30] have been used. The dashed straight line indicates the 
resonance energy of 0.34 MeV. 



of the method adopted for the computations. This does 
of course assume that the methods are mathematically 
and/or physically equivalent which not always is easy to 
confirm. 



A. Interaction dependence 

An asymmetric and numerically more difficult system 
is found in 17 Ne ( 15 0+p+p), which provides an example 
of an exceedingly long-lived two-proton decaying state. 
This nucleus has a resonance with spin and parity 3/2~ 
with an energy of 0.34 MeV above threshold. The ground 
state of 16 F ( 15 0+p) is a resonance at 0.53 MeV above 
threshold. Decay of the 3/2~ resonance in 17 Ne via a 16 F 
resonance is then not allowed due to energy conservation. 
Also, core mass and charge are substantially larger than 
the values of the valence particles, meaning that the es- 
timates of Eqs.(10) are applicable. From them we get 
Pt ~ 96 fm and £ x < 6. From Eq.(16) we then estimate 
that K max « pt/(l fm), i.e. to reach an accurate turning 
point when the two-body channels are exploited we can 
expect K max -va\ues of about 100. 

The numerical difficulties are related both to less sym- 
metry due to the presence of non-identical particles and 
the non-zero intrinsic spins of 1/2 of all three particles. 
The number of partial waves are then immediately rather 
large due to the many allowed spin couplings. The most 
important interactions for the resonance structure are the 
s, p and d-waves of the proton- 15 system, each produc- 
ing two low-lying resonances in 16 F. 



The lowest adiabatic effective potential for the 3/2~ 
resonance is shown in Fig. 4. Except for the last case in 
the figure legend (thick-dotted curve) all the calculations 
have been performed using the interactions in [10]. Wc 
notice the same pattern as in Fig. 3, i.e. increasing K max 
up to values of about 70 is required for convergence of 
the effective potential up to the outer turning point. This 
turning point appears at about 90 fm, in agreement with 
the estimation from Eq.(10). Otherwise the barrier is 
overestimated, giving rise to too small resonance widths. 
Also, the number of partial waves needed to obtain a 
proper convergence is relatively small (£ x and £ y values 
not bigger than 5 is enough, also in agreement with the 
estimation from Eq.(10)). If only s, p, and ci-waves are 
used, the barrier (as illustrated by the thin-dotted curve) 
is too small, and consequently the computed width is 
too big. This barrier is the one used in ref.[10] where a 
resonance width of 3.6 • 10~ 12 MeV was obtained. The 
true width is therefore expected to be a few orders of 
magnitude smaller than this number. 

The thick dotted curve in Fig. 4 has been obtained with 
K m ax = 20 and including all the partial waves consis- 
tent with this value, in total corresponding to the same 
Hilbcrt space as in [30]. This curve is also obtained 
with the same interactions as in [30] where a width of 
4.1 • 10~ 16 MeV is quoted. This width is consistent with 
the 5.4 • 10~ 16 MeV obtained with the WKB approxima- 
tion for this potential. However, as seen in the figure the 
barrier is substantially overestimated compared to the 
results with a much larger basis. The implication is that 
the width should be bigger than this number. 

In Fig. 4, the thick solid line corresponds to a calcu- 
lation with £ x ,£ y < 9, and K max =70 for all the partial 
waves, except for those of large amplitude where K max 
values between 130 and 160 are used. Then the potential 
barrier shows the correct convergence properties and the 
WKB width for this potential and the resonance energy 
of 0.34 MeV is found to be 1.7 • 10~ 14 MeV. This value 
is, as expected, in between the two given in [10] and [30], 
and consistent with the corrected width of (5 — 8) • 10~ 15 
MeV obtained in [29] for a larger K max value. As seen 
in the figure, an increase of K max up to 40 improves sig- 
nificantly the calculated effective potential, although the 
barrier is still a bit overestimated. It it important to 
keep in mind the limitations inherent to the WKB ap- 
proximation. Only one adiabatic potential is included, 
and additional inaccuracy arises from the definition used 
for the knocking rate and the fact that a preformation 
factor of unity is used. 

The potentials shown in Fig. 4 are to a large extent 
determined by the Coulomb and centrifugal potentials. 
This is understandable, since many of the crucial proper- 
ties are determined at distances larger than the ranges of 
the short-range interactions. To investigate the depen- 
dence on these two-body interactions we can compare 
the results arising from the two potentials obtained with 
K max =20. The first of them (long-dashed-dot curve) has 
been obtained with the proton- 15 interaction given in 
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[10], while the second one (thick-dot curve) uses the one 
in [30]. These two interactions are very different, espe- 
cially in their parametrization of the spin-dependence. 
In [10] the spin-spin and spin-orbit operators are s c ■ j p 
and i ■ s p , where j p = £ + s p , s c and s p are the spins 
of the core and the proton, and £ is their relative orbital 
angular momentum. In [30] the more symmetric form of 
s c ■ s p and £ ■ (s c + s p ) is used. 

The second set of spin-dependent operators in [30] docs 
not preserve the usual mean field quantum numbers cor- 
responding to the core, and is as such inconsistent with 
the description of the protons in O. This symmetry 
breaking is especially problematic in cases where one 
spin-orbit partner is occupied by core nucleons while the 
other is available for the valence nucleons. Then the d$/2 
and d 3 / 2 resonances in 17 Ne are mixed in the two-body 
description of the 1~ and 2~ resonances of 16 F. In the 
three-body problem of 17 Ne the two valence protons are 
then forced to partly occupy the same orbits. The Pauli 
principle between core and valence protons is violated 
and unwanted properties may appear, see [40] for details. 
To compensate for these effects a fully phenomenological 
J 71 " dependent three-body potentials are used but their 
effects on the widths could be rather unpredictable. An 
example of unwanted properties is seen in [30] , where the 
potentials obtained with symmetry breaking spin opera- 
tors are adjusted with three-body potentials to reproduce 
the known two-body spectrum in 16 F. However, this si- 
multaneously results in an additional 2~ resonance in 1 F 
at 2.8 MeV with a width of 0.26 MeV. This resonance and 
the corresponding bound state in 16 N at —0.38 MeV are 
not mentioned in the description of the interaction [30] , 
and furthermore they are not known experimentally 

In Fig. 4 it is shown how such exceedingly different two- 
body interactions lead to effective potentials in almost 
perfect agreement, provided that the same basis is used. 
Only a small difference is found around the inner turning 
point which is more sensitive to properties of the short- 
range interactions. This strongly indicates that the spe- 
cific values of the energies of the two-body resonances are 
the only decisive quantities for the intermediate distances 
corresponding to the lowest adiabatic effective potential. 
Small-distance properties essential for spectroscopic fac- 
tors as well as dynamic evolution of the resonance struc- 
tures with distance due to couplings between potentials 
could in contrast be sensitive to the design of the interac- 
tions. These properties are determined inside or outside 
the barrier, respectively. 

The conclusion from these examples is that the width 
is surprisingly insensitive to the two-body interactions as 
long as they provide the proper attraction for the low- 
est adiabatic potentials. This conclusion has two crucial 
assumptions, i.e. first the three-body resonance energy 
is adjusted to the correct energy by use of a short-range 
three-body potential. Second, the effective barrier is ac- 
curately computed for example by use of a sufficiently 
large basis where it is important to allow the higher par- 
tial waves in the Hilbert space although the correspond- 



ing two-body interactions do not have to be precise. The 
explanation is simply that Coulomb potential and cen- 
trifugal barriers dominate at the intermediate distances 
where the confining barrier is located and the width in 
turn is determined. In cases where this turns out to be in- 
correct the width may depend much more on the specific 
choice of the two-body interactions. The most tempting 
guess of such a situation is a large width corresponding 
to a narrow barrier and an outer turning point at a rel- 
atively small distance. 



B. Method dependence 

The methods to calculate resonances defined as poles 
of the S'-matrix must all make use of analytic continua- 
tion into the complex plane. This inevitably involves an 
approximation to the physical quantity of interest, e.g. 
the imaginary part or the resonance width is not found 
in a physical process and therefore not directly an ob- 
servable. The connection has to be established through 
theoretical derivations and model dependent interpreta- 
tions. The continuation is only possible when an ana- 
lytical form is available as for example the potentials in 
the Schrodingcr equation. If they are given as numeri- 
cal tables obtained by fits to measured cross sections at 
discrete points, the numbers must be connected by an 
analytical expression. This is the same result as if an- 
alytical parametrized potentials from the beginning are 
adjusted to reproduce experimental values. Either way 
the potentials can be continued into the complex plane. 

The further away from the real axis, the larger is the 
uncertainty arising both from the model dependent inter- 
pretations in terms of observables and from the somewhat 
arbitrary choice of the initial analytical form. Thus large 
computed widths are for these reasons intrinsically more 
uncertain than the small ones. These methods are equiv- 
alent if the same potentials are employed. The choice 
of method is then only a matter of numerical, and per- 
haps mathematical, convenience. One possibility is the 
use of complex scaled coordinates, i.e. all lengths in 
the Schrodinger (or Faddeev) equation are multiplied by 
cxp(z#), where 9 is a given angle which has to be larger 
than half of the angle of rotation from the real energy 
axis to the direction of the resonance defined in the com- 
plex energy plane. Then the resonance wave function in 
the rotated space is an eigenfunction with bound state 
boundary conditions for the corresponding complex en- 
ergy [26]. 

Rotating the resonance solution back to the real coor- 
dinates results in a wave function with only outgoing flux 
and the same complex energy. This boundary condition 
and a complex energy without complex coordinates are 
then fully equivalent. A third method is to continue the 
interactions analytically by varying a strength parame- 
ter [27, 28]. The results are then obtained as function of 
this parameter and in the end extrapolated back to the 
correct physical value. An example is the computation of 
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the broad + -resonance in 12 C [41]. However, the large 
width is already an indication of an inherent uncertainty. 
For all these methods the computed widths are model 
dependent approximations most accurate for relatively 
small widths. 

An almost identical method is the Gamow shell model 
which so far only is applicable for nucleons and a core 
[21, 42]. It is the complex rotated mean-field shell model 
with wave functions expanded on an ordinary single- 
particle basis. The computed eigenvalues are investi- 
gated as functions of the number of basis states and their 
spatial extension. A few of the eigenvalues converge to 
specific complex energies corresponding to resonance po- 
sitions and widths, while the majority only reflect an 
attempt to discretize the continuum. Only uncorrelated 
motion can be described. 

It is also possible to remain on the real axis, both 
energy and coordinates, but the resonance energy and 
width must be defined in a different way. In a scattering 
formulation the elastic cross section varies through a peak 
over a small range of energy. Apart from quantitatively 
important and qualitatively unimportant corrections, the 
peak position and its width are the resonance energy and 
the width. A large width then means that the cross sec- 
tion is smeared out and resembling the background which 
obviously makes both energy and width determination 
more uncertain. 

Another definition of the width is through the time de- 
pendent, or stationary wave packet of an evolving wave 
function. An initial condition then has to be assumed, 
e.g. a source term in the Schrodinger equation adding 
probability at small distances. The decay rate by which 
this probability disappears is then directly interpreted as 
the resonance width. When more resonances contribute 
this is much more complicated and turns into a coupled 
channels problem. Analyses of experimental data use this 
formulation expressed as transition probabilities through 
potential barriers corresponding to two consecutive co- 
herent two-body decays [43]. 

The simplest version of a stationary incident wave 
packet attempting to tunnel through a barrier is esti- 
mated by the semiclassical WKB-tunncling probability. 
This formulation is still possible for a multi-channel prob- 
lem, simply by choosing a one-dimensional path through 
the many-dimensional space. By definition this is a semi- 
classical method where the tunneling probability should 
be small and the potentials should be smooth. 



C. Perturbation treatment 

Another possibility is the perturbation treatment used 
in [11]. The width is obtained in three steps. First the 
bound state problem is solved in a box of hyperradius p 
less than the outer turning point p t . Second, this wave 
function is used as a source term of arbitrary strength 
to find the resonance wave function, and third the decay 
rate (or the width) is found by computing the outgoing 



flux at an arbitrary large p max distance. At this stage it 
is then necessary to impose the proper boundary condi- 
tion to the computed wave function, which for a system 
of three charged particles is not known explicitly A de- 
tailed discussion about possible different approximations 
to implement the Coulomb boundary condition in the hy- 
perspherical harmonic expansion method can be found in 
[44] . In our adiabatic approach it is enough to impose to 
the radial solution to go as G^(i], p) + iT^T], p), where T 
and Q are the Coulomb functions, and the index £ and the 
Sommerfeld parameter r\ are obtained numerically from 
the adiabatic potential (see [45] for details). 

The assumption is that the perturbatively computed 
wave function describes the true resonance, and the re- 
lated outgoing flux at the corresponding p max radius 
gives the decay rate. The first of the steps does not re- 
quire knowledge of the details of the outer part of the 
barrier. However, calculation of the outgoing flux at a 
given p m ax value beyond the barrier obviously requires 
barrier knowledge at least up to p ma x- 

Therefore, for this procedure to work it is necessary 
either to compute accurately the potential barrier up to 
Pmax or to make some assumption about the neglected 
outer parts of the barrier, which still are inside the turn- 
ing point. This is unavoidable since an infinitely thick 
barrier smoothly continued from the box radius p inside 
the turning point would give vanishing decay rate (width 
equal to zero). The opposite assumption of vanishing 
barrier outside the box would result in a finite decay rate. 
This means that if the barrier is overestimated, as for too 
small K max , the outg will be reduced 

producing a too small width. It is then meaningless to 
expect an accurate computation of the tunneling prob- 
ability through a barrier without accurate knowledge of 
that barrier. 

In the same way, one should be very careful when com- 
puting observables obtained from the behaviour of the 
wave function at very large values of the hyperradius (as 
for instance 1000 fm, as quoted in [29, 30]). Direct com- 
putation of these large- distance properties is in obvious 
conflict with the poor knowledge of the outer barrier. 
If typically one needs a large value of K max to get ac- 
curate calculations up to the outer turning point, this 
K m ax value should be clearly bigger if accurate results 
are required for distances of a few times p t - 

The procedure applied in [8, 9] circumvent this prob- 
lem of direct calculation at unreachable large distances. 
First it is tempting to employ momentum space instead of 
coordinate space but this would at best only change the 
problem into uncertainties at large momenta or equiva- 
lently the small distances decisive for the resonance struc- 
ture would be uncertain. If both small and large dis- 
tances are needed accurately it is equally convenient to 
work in coordinate space where the small distances nat- 
urally are most accurately computed. To get accurate 
large-distance asymptotic properties, the desired observ- 
ables should be computed with acceptable accuracy for 
a given large basis at the largest possible distance. The 
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computed observable must then remain unchanged when 
the distance is further increased provided the basis is 
correspondingly increased. Passing this test is equivalent 
to sufficient accuracy at both small and asymptotically 
large distances. 

Comparison of different methods is quantified in ta- 
ble I where the third column, T pert , gives the widths of 
the Hoyle resonance in 12 C computed using the pertur- 
bative method described in [11] with inclusion of only 
the first adiabatic potential, precisely as for the WKB 
results, Twkb, in the second column. We note that the 
perturbative results are consistent with the ones obtained 
in the WKB approach, although they tend to be about 
one order of magnitude smaller, especially when K max is 
far from the one required for convergence. Again, calcu- 
lations made with a too small value of K max give rise to 
too small widths, while the calculation made with a suffi- 
ciently large K max and a sufficiently large number of par- 
tial waves provides a result in good agreement with the 
WKB estimate and the experimental value. The values 
obtained with the perturbative method depend weakly 
on the p m ax value and the value of V used to obtain the 
wave function with outgoing boundary condition. In any 
case, the computed results do not change significantly. 

Similar agreement is obtained for the 3/2~ state in 
17 Ne. In this case the values of T pert oscillate between the 
4.9 • 1CT 17 MeV obtained with K max =20 (long-dashed- 
dot curve in Fig. 4), and the 3.0 • 10~ 15 MeV obtained 
with K max < 70 (thick solid line in the figure). For these 
two cases the WKB estimates were 5.4 • 10~ 16 MeV and 
1.7 • 10~ 14 MeV, respectively. 

Therefore, the different methods give rise to similar 
widths, at least within the method uncertainties. Once 
the three-body energy of the resonance is given correctly 
(with the help for instance of an effective three-body 
force), the two methods are equally sensitive to the defi- 
ciencies of the potential barriers arising from the use of 
a poor basis. Use of a too small K max value or too few 
partial waves produces incorrect barriers, and therefore 
uncontrolled widths. 



D. Fragment energy distribution 

After three-body decay of a given resonance, one of 
the most investigated observables is the energy distribu- 
tions of the fragments after the decay. As shown in [46], 
these quantities are mainly determined from the prop- 
erties of the coordinate space wave functions at large 
distances. This happens because the hyperspherical har- 
monics transform into themselves after Fourier transfor- 
mation into momentum space. The kinetic energy distri- 
bution of the fragments at a given p is therefore, except 
for a phase-space factor, obtained as the absolute square 
of the total wave function in coordinate space for that 
p, but where the five hyperangles are interpreted as in 
momentum space. 

To obtain the energy distributions after decay one ob- 




FIG. 5: Proton energy distributions for the 3/2~ resonance 
in 17 Ne where the maximum allowed emission energy is used 
as unit. The left part (a) is for p=25 fm for different values of 
K max - The right part (b) is for a very large basis for increas- 
ing hyperradii. Only the lowest adiabatic effective potential 
is used. 



viously needs to consider a value of p clearly beyond the 
potential barrier, which means values of the hyperradius 
of several times the outer turning point. Following the 
discussion in the previous section, it is clear then that 
an accurate computation of such distributions requires a 
basis able to describe properly the three-body state at 
such large distance. In other words, reliable calculations 
require sufficiently large values of K max , which will be 
larger than the ones reproducing the potential barrier up 
to the outer turning point. 

For a system like 17 Ne, where the outer turning point 
for the 3/2" resonance is about 90 fm, the basis needed 
to compute the energy distributions after decay becomes 
soon very big, and therefore the numerical calculation 
very heavy. However, to illustrate the importance of a 
sufficient size of the basis, we first show in left part of 
Fig. 5 the energy distribution of the proton when differ- 
ent K max values are used and when p is fixed at 25 fm. 
For such p value, as seen in Fig. 4, accurate calculations 
are possible. Although this value of p is clearly below 
the outer turning point, it illustrates the dependence of 
the energy distributions on the basis size. The only thing 
to have in mind is that the plotted distributions corre- 
spond to an intermediate stage, before the decay of the 
resonance. As a consequence, the plotted distributions 
do not necessarily resemble the experimental ones, which 
should be computed at a much larger value of p. The 
interplay between the two-body interactions can clearly 
change the energy distributions when moving from small 
to large values of p. 

In any case, from Fig. 4, we know that at p=25 fm a 
K max value of at least 40 is needed to have a converged 
potential at this distance. This fact is also reflected in left 
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part of Fig. 5, where we see that the energy distribution 
obtained when K max =40 matches pretty well the result 
obtained with a much bigger basis (thick solid line) . This 
energy distribution corresponds to a situation in which 
one of the protons takes most of the energy, while the 
other one stays close to the core. However, when the ba- 
sis size decreases the proton energy distribution begins 
to fill the intermediate region, such that when K max —20 
(thin solid line) only one wide peak at about 0.5 is seen. 
The same strong dependence on K max as the one ob- 
served in the figure can be found for larger values of p, 
in particular for the values needed when computing the 
energy distributions after the decay. 

To see the dependence on the hyperradius we show in 
the right part of Fig. 5 how the proton energy distribution 
changes with the largest basis used in the left part. Now 
increasing p shifts the proton energy distribution from 
two narrow peaks at low and high energy to one broad 
peak around half the maximum energy. These distribu- 
tions have not yet converged as function of p. A larger 
value is needed with a correspondingly larger basis but 
the trend is clear. The protons are emitted roughly with 
equal energy. It is very illuminating to notice that the 
correctly converged result obtained for a large hyperra- 
dius and a large basis resembles an inaccurate result from 
a small basis and a much too small hyperradius. 

VI. SUMMARY AND CONCLUSIONS 

The partial width for a resonance decaying into three 
fragments must necessarily deal with the corresponding 
three-body problem. As for a-emission we believe the 
crucial ingredients are provided by effective potentials at 
intermediate and large distances measured relative to the 
radius of the decaying nucleus. The many-body effects 
essentially only enter at small distances as preformation 
or spectroscopic factors and in microscopic derivations of 
the effective two-body interactions. 

We concentrate in this paper on the three-body prob- 
lem. We assume that the small distance boundary condi- 
tions necessary for resonance computations is simulated 
by the correct three-body energy and an attractive pocket 
within the confining barrier. We furthermore assume it 
is possible to specify the effective two and three-body in- 
teractions. With these assumptions we have isolated the 
crucial three-body problem from the underlying many- 
body degrees of freedom. Accurate resonance structures, 
partial decay widths, and fragment momentum distribu- 
tions must first of all be computed for this well defined 



few-body problem. Inclusion of many-body effects may 
be important but this cannot avoid the three-body prob- 
lem. The more elaborate schemes only add to the com- 
plexity involved in obtaining three-body input parame- 
ters. 

We first formulate the basic concept of effective three- 
body potentials. As a test we then compare experimental 
widths for a number of 12 C resonances with results from 
the simplest WKB application. As for a-emission the 
largest variation is reproduced only leaving effects from 
preformation factors or equivalently many-body effects 
amounting to less than one to two orders of magnitude. 
The reproduced systematics is not monotonous with ex- 
citation energy or angular momentum. 

The essence of the problem is now narrowed down to 
accurate computations of the effective three-body poten- 
tials. The classical turning points for the dominating 
potential therefore specify the region of interest in coor- 
dinate space. Any uncertainty in this region is enhanced 
exponentially in the computed widths. We give analyti- 
cal estimates of both turning points and the number of 
contributing partial waves. We also estimate analytically, 
and demonstrate numerically, which basis size is neces- 
sary for accurate computations of the potentials in this 
coordinate region. Almost all published results employ 
insufficient basis sets. 

We investigate the dependence of the computed widths 
on the effective two-body interactions. As for a-emission 
the total energy is crucial and the two-body interactions 
should only provide roughly the same resonance energies 
independent of contributing spin and orbital angular mo- 
mentum structure. The decisive properties are supplied 
by Coulomb and centrifugal forces. We discuss the differ- 
ent methods applied to width computations and conclude 
that they in most cases are equivalent. 

Finally we illustrate how it is more difficult to get 
accurate momentum distributions of the fragments af- 
ter three-body decay. These observables are sensitive 
to properties of the coordinate space wave functions at 
distances outside the turning point. Therefore an even 
larger basis is required. Unfortunately, a two-fold in- 
accurate computation with a too small basis and a too 
small hyperradius by a strange coincidence resembles the 
correctly converged result. 

In conclusion, the partial three-body widths for decay 
of many-body resonances are first of all determined by 
three-body properties. Many-body effects are less impor- 
tant. The effective three-body potentials or equivalently 
the corresponding three-body wave functions are decisive 
and must be accurately computed. 
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